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We present preliminary results for the mass of the 6q flavor-singlet state (J^ =0^, S = —2) called If-dibaryon, 
calculated in quenched QCD on 16^x30 and 24^x30 lattices with improved gauge and fermion actions (Symanzik 
improvement, Clover action). For both lattice sizes we applied the fuzzing technique to enhance the overlap 
with the ground state. We observe a H-mass above the AA-threshold for strong decay. The difference in mass, 
niH — 2mA, increases with increasing lattice size. 



1. INTRODUCTION 

In 1977 a bound six-quark state {uuddss), the 
H-dibaryon, was predicted in a bag-model calcu- 
lation by Jaffe [1]. This state is the lowest SU(3) 
flavor singlet state with spin zero, strangeness -2 
and = 0+. These dibaryons may play an im- 
portant role as Bose condensate in nuclear mat- 
ter prior to the quark-hadron phase transition of 
QCD at high density. 

In the last twenty years many attempts to verify 
the existence and stability of this particle were 
undertaken by means of various methods. Per- 
turbative calculations included spin-dependent 
q-q interactions arising from one gluon exchange 
(OGE) [1,2], instanton induced interactions (III) 
[3,4] and Goldstone boson exchange (GBE) [5]. 
The attractive q-q interaction is controlled via 
the potential Vs ^ (A°A2)(siS2), where A° denote 
SU(3) color (OGE) or flavor (III/GBE) genera- 
tors. The resulting H-masses scattered in a range 
of a few hundred MeV around the 2231 MeV 
AA-threshold for strong decay. 

The lattice attempts to calculate rriH a decade 
ago [6,7] gave contradicting results, but a new 
calculation [8] showed that the q-q attraction 
decreases with increasing volume. Hence it has 
been suggested that the H-dibaryon is unbound 
(approx. 110 MeV above the threshold) in the 
infinite volume extrapolation. 
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2. DETAILS OF THE SIMULATION 

Our spectrum calculation was performed in 
quenched QCD with Wilson fermions on lattices 
of the size 16'^ x 30 and 24^ x 30. We gener- 
ated 57 independent gauge field configurations 
for the smaller lattice and 15 configurations so 
far for the larger lattice, separated by 100 sweeps 
of 4 overrelaxation and one heatbath step each. 
For the gauge sector we applied the tree-level im- 
proved (1,2) Symanzik action at a gauge coupling 
P = 4.1, which corresponds to a lattice spacing 
a — 0.177(8) fm defined by the string tension [9]. 
A tree-level improved clover term was used in the 
fermionic part of the action. 

The correlation functions were calculated via the 
relation C{t) = {O{0)OHt)) with appropriate 
operators for the lambda and the H-dibaryon 
[10] (color indices: roman letters, spinor indices: 
greek letters): 

Oa{x) = eabc{Cj5)04Kix)d'l3{x)Syix) 

+d^^ (x)4 (x)u^; (x)- 2sl (x)4 (x)d; {x)] 
Oh{x) ~ 3{udsuds) — 3{ussudd) — 3{dssduu) 

(abcdef) = eabc^defiC-f5)apiCj5)^siCj5)e<P 

*a^(x)b^(a:)<(x)d^(x)eKa:)f^(x) 

We summed over spatial coordinates to project 
onto zero momentum and over all spinor and color 
indices to get diagonal correlators. 
Quark propagators were calculated at five differ- 
ent values of the bare quark mass, which have 
been estimated to be in the range 50 - 250 MeV. 
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Figure 1. Influence of fuzzing on the dibaryon 
correlation function through different number of 
fuzzed quark propagators 

The mass is controlled through the hopping pa- 
rameter K„ for the degenerated u- and d-quark 
mass and Kg for the s-quark mass. The physical 
K„ used in the extrapolation was set by the ratio 
of nucleon mass and string tension. The value for 
Kg was determined by the mass ratio of lambda 
and nucleon. 

In addition to the H-dibaryon and lambda we cal- 
culated the correlation functions for the strange 
mesons K, K* and the S baryon. The obtained 
particle masses lie all in a 10% range around the 
experimental values, which is a common result in 
calculations using the quenched approximation. 

3. FUZZING 

Extended operators and smeared gauge fields 
provide a better overlap with the ground state 
[11]. For the quark fields we use: 

R) U\x - ij)...U\x - RfL)^{x - Rfi) 

+ U{x)...U{x + {R- l)/i)*(.T + Rfl)) 

Smeared gauge fields improve the simulation of 
the gluon cloud around the quarks: 

U'{x,fi) ^VsU{3)[cU{x,fl) + ^ Ustaples(,^)] 

The application of fuzzing for two of the six 
quarks inside the dibaryon flattens the curvature 
of the effective mass (Fig. 1). The largest plateau 
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Figure 2. Lambda and H-masses extrapolated to 
the physical 0.149 (line indicates physical Kg) 

in the region with small errors is obtained with 
fuzzed quarks for the four lighter u- and d-quarks. 
Therefore we used this variant to calculate our 
correlation functions. 

4. RESULTS 

4.1. Comparison of Dibaryon and Lambda 

Masses 

The particle masses were extracted from the 
long range behavior of the correlation functions. 
We leave out successively data points in the ex- 
ponential fits at small t to yield stable results. 
Contributions of excited states were substantially 
reduced by the fuzzing technique. Additionally, 
fits with two exponentials enlarge the plateau 
even for the first time slices. 

The results on the smaller lattice show a H-mass 
in the same range as twice the A-mass (Fig. 2). 
We obtained tuh = 2221(141) MeV, which is just 
below the AA-threshold for strong decay. For 
the difference in mass we obtain rriH — 2m\ = 
— 10 ± 141 MeV. It is therefore impossible to de- 
cide, whether the dibaryon is unbound or slightly 
bound. Moreover, we cannot rule out consid- 
erable finite size effects on the small (2.8 fm)"^ 
lattice. 

The mass splitting on the larger lattice (4.2 fm)'^ 
seems to be more pronounced (Fig. 3). The H- 
mass is larger than two A's for the parameter 
combinations calculated so far. Although our 
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Figure 3. Lambda and dibaryon masses on both 
lattices for k„ = 0.1475 

statistics on the large lattice is not yet sufficient 
to draw a definite conclusion, the preliminary re- 
sults shown in Fig. 3 suggest that the H-dibaryon 
is unbound in the infinite volume limit. 

4.2. Ratio of the Correlation Functions 

The errors in the correlation functions for the 
dibaryon and the lambda are strongly correlated. 
Some of the statistical uncertainties can thus be 
eliminated by analysing directly the ratio of corre- 
lation functions. The ratio of the correlation func- 
tions of dibaryon and lambda should rise, when 
niH < 2toa: 

. C{H) 



-mnt 



C{H).-. _ ^^^^^ . 

or should be constant, when m/f > 2mA: 

C(77)^e-^-* =^ const. 

On the smaller lattice this ratio is rising at i > 8 
with and without fuzzing (Fig. 4), which pro- 
vides evidence of a slightly bound H-dibaryon. 
On the larger lattice the ratio remains constant 
until i = 8. With our present statistics the sig- 
nal gets lost for larger values of t. Nonetheless, 
the large plateau of constant values again sug- 
gests that we see an unbound H-dibaryon on this 
larger lattice. 

5. CONCLUSIONS & OUTLOOK 

With our present, still preliminary, analysis of 
the H-dibaryon and A masses we confirm the find- 



Figure 4. Ratio of the correlation functions for 
K„ = 0.147 (0.1475 for 24^x30) and Ks = 0.146 

ings of Negele et al. [8] presented at last year's lat- 
tice conference. We observe that the difference in 
mass, vtlh — '^vnt^^ increases with increasing lat- 
tice size and therefore will lead to an unbound 
H-dibaryon in the infinite volume limit. 
We will continue our simulations on the larger 
lattice with more configurations and additional 
combinations of k„ and Ks to extrapolate to the 
physical values of the masses. Simulations on a 
smaller lattice would be advantageous to confirm 
the connection between finite lattice size and the 
size of the mass splitting mu — 2m\. 
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